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Abstract 

We suggest that low-lying eigenvalues of realistic quantum many-body hamil- 
tonians, given, as in the nuclear shell model, by large matrices, can be calcu- 
lated, instead of the full diagonalization, by the diagonalization of small trun- 
cated matrices with the exponential extrapolation of the results. We show 
numerical data confirming this conjecture. We argue that the exponential 
convergence in an appropriate basis may be a generic feature of complicated 
("chaotic") systems where the wave functions are localized in this basis. 
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Statistical properties of complex quantum systems have been studied extensively from 
various viewpoints. The seminal papers by Wigner and Dyson developed the random 
matrix theory (RMT 0]) where the systems are considered as members of a statistical ensem- 
ble, and all hamiltonians of given global symmetry appear with certain probabilities. The 
canonical Gaussian ensembles correspond to systems with very complicated dynamics 
when, in almost all bases connected by the transformations preserving global symmetry, 
the components of generic eigenfunctions are uniformly distributed on the unit sphere in 
multidimensional Hilbert space. On local scale, Gaussian ensembles predict specific corre- 
lations and fluctuations of spectral properties which are in agreement with empirical data 
for atoms, nuclei 0, quantum dots and resonators (microwave P] and acoustic ex- 
periments). These spectral features are considered usually as signatures of quantum chaos 



Recently, the detail studies of highly excited states in realistic atomic |]T2[ and nuclear 
||r3| calculations demonstrated that such many-body systems are close to the RMT limit 
although they reveal some deviations, partly due to the presence of the mean field 



coherent components |T5[ of the residual interaction and its two-body character 0,[TB|. In 
complex atoms and nuclei, precise experimental information exists, as a rule, about low- 
lying states only. Effective residual interactions , such as the Wildenthal-Brown interaction 
||T7| for the S(i-shell model turned out to be successful well beyond the input used for their 
original fit. This justifies the use of such interactions for studying generic complicated states 
in the region of high level density. The whole shell model approach is based on the large- 
scale diagonalization even if one is interested in the low-lying states only. The dimensions 
of matrices increase dramatically with the number of valence nucleons which makes the 
full diagonalization impractical, even after projecting out correct angular momentum and 
isospin states. This problem is avoided in the Monte Carlo shell model method WM, but. 



apart from the so-called sign problem [|1^] which requires the introduction of an extrapolation 
when working with realistic interactions, this method is better suited for calculating thermal 
properties or strength functions than spectroscopic characteristics. In order to consider 



individual levels, one needs to supplement the Monte Carlo sampling with some variational 
procedure including an additional "stochastic" diagonalization ( or the QMCD approach 
||2T|). An important step towards larger dimensions in the standard shell model is made with 
the development of the DUPSM code . 



Here we suggest a simple approach which allows one, for calculating energies of low lying 
states, to reduce large dimensions of matrices under study by orders of magnitude, while 
keeping high precision of the results. The approach is based on the statistical properties 



of complicated many-body states |l^,|T3|. Because of the strong residual interaction and 
"geometric chaoticity" ||l3l of the angular momentum coupling, the eigenstates are extremely 
complex superpositions of independent particle Slater determinants. However, in contrast to 
the limiting case of the Gaussian orthogonal ensemble (GOE), the stationary wave functions 
are not fully delocalized in shell model space. Due to inherently self-consistent nature of the 
residual interaction (even if it is extracted in a semiempirical manner), its strength does not 
exceed the typical spacings between single-particle levels which are determined by the mean 
field, i.e. by the same original forces. Together with the fact that the two-body forces cannot 
couple very distinct configurations, this leads to a band-like structure of the hamiltonian 
matrices in the shell model basis. 

Theory of banded random matrices did not reach the same degree of completeness as 
that of canonical Gaussian ensembles. Nevertheless, both mathematical pS) and numerical 



|2J] arguments favor the localization of the eigenstates in Hilbert space, similar to the 
coordinate localization of electronic states in disordered solids. The generic many-body 
states in complex atoms or nuclei have a typical localization width [0,^,^. Inversely, the 



simple shell model configurations are packets of the eigenstates. Their strength function is 
fragmented over the range of energies characterized by the spreading width F which is nearly 
constant along the spectrum because the coupling matrix elements between the complicated 
states are small just as it is needed to compensate small level spacings in the region of high 
level density pB|-p5|. The qualitative arguments are confirmed by more general theory |T3 



as well as by detailed numerical calculations for atoms and nuclei ||T3|j29[| . The nuclear 



case is close to the strong coupling limit [^,^,0 where the typical width can be estimated 
p9| as r ^ 2a in terms of the energy dispersion of a simple configuration \k), 

al = {k\{H-{k\H\k)r\k) = Y.\HM\'. (1) 

Here Hi^i are the off-diagonal elements of the residual interaction between the basis states 
so that the calculation of (P does not require any diagonalization. The dispersions ak of 
different simple states fluctuate weakly |jl3| and in our estimate of T they are substituted by 



the appropriate mean value a which can be found by the methods of statistical spectroscopy 



The practical method of truncating large shell model matrices was suggested in . 
The shell model states are grouped into partitions (sets of states belonging to the same 
particle configuration). Since the states separated in energy by an interval broader than F 
are not significantly mixed with the studied state, we truncated the matrix retaining only 



the partitions whose statistical centroids E = {k\H\k) are closer than 3cr. The spin-isospin 
projection and the elimination of the center-of-mass admixtures can be done within the 
truncated subspace only. In order to keep the correct shell model structure, the partitions 
should be included as a whole. As shown in [Q, this method allows for the calculation 
of low-lying energies with sufficient precision in large shell-model spaces. The truncated 
eigenvectors overlap with the exact ones on the level of better than 90%. 

Going beyond the simple truncation, we consider the convergence of level energies to 
the exact values as a function of the increasing dimension n of the diagonalized matrix. 
As an example we take the ^^Sc nucleus where the p/-shell model dimensions of l/2~ and 
3/2^ states are 13016 and 24474, respectively. Spectroscopic information on this radioactive 
isotope is only tentative providing an interesting experimental and theoretical problem. For 



the calculations, the FPD6 interaction |3^ was used. Fig. 1 shows the calculated energies 
of the two lowest 3/2~ states and two lowest 1/2^ states for several values of n ranging 
from n = 2000 to the full dimension A^. Already the smallest dimensions lead to a good 
agreement within few As the dimension increases, in all four cases the running eigenvalue 



converges very fast and monotonously to the exact result. The convergence is almost pure 
exponential, E{n) = + ylexp(— 7^); typically A SOOkeV, 7 ^ 

The exponential convergence of eigenvalues would be extremely helpful for shell model 
practitioners. It would make almost redundant the full large scale diagonalization if one is 
interested in the low-lying states only. Instead, the calculations for several increasing dimen- 
sions (still far from the full value and therefore easily tractable) would end in determining 
the exponential parameters and simple extrapolation to the exact result. At present, the 
rigorous mathematical theory of convergence is absent, and we limit ourselves by qualitative 
arguments and plausible conjectures. 

The convergence under consecutive truncations is determined by the type of the matrix 
and by the original unperturbed basis which orders the basis vectors in a certain way. 
The ordering is done almost uniquely in the spherical shell model where the mean field 
is fixed and all many-body states are organized in partitions. For the lowest levels, the 
admixtures of highly excited states outside of the starting truncation correspond to the 
wings of the strength function. As confirmed by the atomic and nuclear |r5| , |25|] studies. 



the strength function has in average universal shape. This shape evolves from the standard 
Breit-Wigner function |^ for the "weak damping" case to the Gaussian form at strong 



damping (semicircle in the RMT limit of the uniform spectra [^,^ which is not reached 
in practical cases). Correspondingly, the dependence of the spreading width on the strength 
of the residual interaction changes from quadratic in the standard golden rule [^] to linear 



31| , p!5|p9| . The remote wings of the strength function have their own energy behavior p!2| , p9 



With high accuracy they can be described p9| by an exponential function of energy. This 



is a clear manifestation of the localization of the eigenf unctions typical for the banded 
hamiltonian structures p6| , |23[] . In this limiting regime, the average local strength is simply 



proportional to the total remaining strength, F[E) ^ const [1 — /'^ dE' F[E')] which gives 
F{E) ~ exp(— const E'). The exponentially weak mixing should lead to the exponentially 
small energy shifts and to the corresponding convergence of the eigenvalues. 

We can expect the exponential behavior to be generic for the large matrices of quasi- 



banded form with the off-diagonal elements of approximately the same order of magnitude 
along the spectrum, This conjecture can be checked by generating random matrices with 
the desired properties and diagonalizing them in a sequence of progressing truncations (the 
matrices are first ordered according to their diagonal elements). Since there is no "vertical" 
structure in such random matrices, we do not have here physical arguments concerning the 
optimal truncation sequence. We show in Fig. 2, left, a typical result for the banded GOE 
matrix with the width b = 0.293A^ which clearly demonstrates the exponential convergence. 
The full GOE matrix. Fig. 2, right, converges more slowly in the absolute sense and does 
not saturate. This is a simple consequence of the fact that the ground state of a GOE matrix 
is repelled by the higher states to the edge of the semicircle (-2 with the GOE definition 



accepted here and in |T3|)- This process is driven by the off-diagonal elements; all of them 
in average have the same order of magnitude. Their number and, whence, the dispersion a, 
eq. (|l|), is greater in the full GOE. Therefore the distance from the unperturbed position 
is also greater in this case (the diagonal and off-diagonal contributions add in quadratures). 
In all studied cases with the width b changing from O.IA^ to the full GOE, the convergence 
is exponential and the exponent 7 is approximately scaled inversely proportional to b. 

In realistic cases there is also a leading sequence of regular diagonal elements. A similar 



banded matrix example, considered in |12[, goes back to Wigner The matrix consists 
of the equidistant diagonal with the spacing D and random off-diagonal matrix elements 
Vki within the band \k — l\ < b. At relatively weak interaction, g = {V'^)/D'^ < 1, the 
main contribution in the perturbation series for the admixture w„ = of a very remote 
state |n), n ^ 1, to the wave function of a low-lying state |0) is given by the summation 
of long "straight" paths in Hilbert space connecting the initial state with the final one 
through various intermediate stops. Because of the random character of the off-diagonal 
interaction, the mean value of Wn is determined by the squares of the contributions of these 
paths (no interference). In the approximation of a weakly changing level density, this can 
be approximately written as an integral equation of the random-walk type. 
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= ^/n + ^ X! fn-kWk (2) 

where the factor comes from the energy denominators, and /„ shows the behav- 

ior of typical squared off-diagonal matrix elements V^q as a function of the distance n 
from the diagonal. With the sharp band boundary ||12[, the weights Wn decrease very 
fast, essentially as exp(— nlnn) ~ (n!)^^. With the smooth cut-off, the convergence is 
getting closer to exponential. Thus, for the exponential cut-off of the matrix elements, 
Vki ~ exp(— |fc — we have /„ = exp(— 2n/6), and eq. (0) allows a simple solution 

Wn = Aexp(— 2n/6)/'n,^. Therefore the contributions to energy of the state |0) should con- 
verge ~ nDwn = DAexp{—2n/b)/n. Fig. 3 illustrates this consideration by an example 
of the numerical diagonalization of a random matrix with the equidistant diagonal and the 
exponential cut-off. One may note that the rate of convergence is similar to that in the 
shell model calculation, see above, where the effective width of the band is close to 6 ~ N/4 
|[T3| . However, the method of eq. (0) becomes invalid in the case of strong interaction when 
the contributions to the perturbation series of additional loops in Hilbert space cannot be 
neglected. The range of convergence is seen from the expression for the constant in the 
above solution, A = g/[l — gJ2k=i k^"^], which determines the critical value Qc = G/vr^. 

It is interesting to test the character of the convergence in simple solvable models. A 
harmonic oscillator, shifted from the equilibrium position by a constant force, H = a^a + 
A (a -|- a'^), lowers its energy by A^. The exact ground state is a coherent combination of 
unperturbed states \n). In accordance with the composition of the coherent state, the 
convergence of the energies in the unperturbed basis of the original oscillator is ~ A^"/n!. 
This is clearly seen in Fig. 4, left. The fast convergence is due to the constant level density 
along the main diagonal while the perturbation has matrix elements growing only ~ ra^/^. 
In the case of a quartic anharmonic oscillator, the exponential convergence is modulated 
with oscillations. 

Another example displays the case of the slow convergence. The tight-binding model 
of a finite one- dimensional lattice has degenerate levels in each of N identical wells and 
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the amplitude v of hopping between the adjacent wells. The eigenstates of the model are 
delocalized standing Bloch waves with energies within the band, Eg = 2vcos(pq, ipg = 
nq/ (A^+ 1), q = 1,2, . . . , N. As it is easy to see, the truncation in the site basis corresponds 
to the convergence ~ l/n^, see Fig. 4, right. 

Tridiagonal matrices with the entries H„n = e„, and Hn-i.n = Hn.n-i = smoothly 
depending on n can be analyzed in a general way using the recurrence relation for the 
secular determinants Dn{E) of the matrix {H — E) truncated at the nth step, 

Dn{E) = (e„ - E)Dn-i{E) - V^Dn-2{E). (3) 

The consecutive approximations E^'^~^^ and E'^'^\ n S> 1, to a low- lying eigenvalue E are 
the roots of = and = 0, respectively. For e„ > E^") (with a slight 

modification, the method works also for initially degenerate matrices with e„ =const), the 
asymptotic behavior of the energy increments A„ = i?*^") — E'^""^) follows from as 

''^-^K. (4) 



(A„ + A„_i)(A„_i + A„_2) e„e„_i 

The exponential convergence corresponds to A„ ^ A = constT^ at large n (similar increase 
of diagonal and off-diagonal elements). Then the increment ratio ^„ = A„/A„_i also goes 
to a constant limit ^ = (1/2A)[1 — (1 — 4A^)^/^] which restricts the exponential convergence 
region to A^ < 1/4. The existence of the constant limit is still compatible with an additional 
preexponential factor weakly dependent on n; at large n corresponding fits are usually 
indistinguishable. An explicitly solvable (by the Bogoliubov transformation) model of the 
harmonic oscillator with the perturbation A(a^ + a^^) agrees completely with this estimate. 
The case A^ = 1/4 corresponds here to the degeneracy of the oscillator with zero frequency, 
and at A^ > 1/4 the spectrum is inverted. In general, A = 1/4 resembles a critical point; 
the convergence here is described by a power law being exponential outside of this region. 

The main difference of the tight-binding case from the oscillator model is the degeneracy 
of the unperturbed levels (absence of the leading diagonal) which results in the delocalized 
wave functions of eigenstates. Because of the degeneracy, the analog of eq. (H) contains, 
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instead of A^, the ratio v'^ / E^"^^ E^"^ which is just equal to the critical value 1/4 in the limit 
of large n. The situation is similar for the spin chains with the nearest neighbor interaction 



where the finite size effects on the ground state energy were repeatedly studied fS^] and 
corrections also go as ra"^. We expect the presence of disorder (random positions of the 
original site levels), which leads to the localization of the eigenfunctions, to be accompanied 
by the transition to the exponential convergence of the eigenvalues. 

In conclusion, we discussed the convergence of the low-lying eigenvalues of large matrices 
describing the realistic many-body hamiltonians of shell-model type. We gave arguments in 
favor of the conjecture that the exact diagonalization of relatively small matrices, truncated 
on the grounds of physical partitions and generic spreading widths of simple configurations, 
provides a starting approximation which can be extrapolated to the exact result with the 
aid of a simple exponential continuation. The arguments are based on the generic features of 
quantum chaotic many-body dynamics, simple models and the results of numerical analysis. 

The authors acknowledge support from the NSF grants 95-12831 and 96-05207. They thank 
F.M.Izrailev for stimulating discussions. 
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FIGURE CAPTIONS 



Figure 1. Energy deviations from the exact shell model result for the lowest excited states 
l/2~ and 3/2~ in ^^Sc (diamonds) calculated as a function of the progressive matrix trunca- 
tion n. Solid lines give a fit Aexp{—jn) with the parameter values A = 0.26, 0.19, 0.23, 0.35 
MeV and 7 = 49.7 x 10-^47 x 10-^25 x 10-^30 x 10"^ for 1/2^, 1/22 , 3/2r, and 3/22, 
respectively. 

Figure 2. Energy deviations for the ground state of random matrices of dimension 
1000 as a function of the progressive matrix truncation n (diamonds): the GOE-hke banded 
matrix of the width b = 0.293N, chosen in such a way that a half of the matrix elements 
vanish, approximately in the same proportion as in typical shell model cases, left panel; 
the full GOE matrix, right panel. Sohd hues show a fit Aexp{—^ri) with A — 1.18 and 
7 = 3.8 X 10-^ left, and A = 1.96, 7 = 1.53 x 10"^ right. Note the absence of the 
horizontal asymptotics in the case of the full matrix. 

Figure 3. Convergence of the ground state energy for the banded random matrix with 
the exponential cut-off of matrix elements Vki ~ exp{—\k—l\/b) with b — 0.293iV (diamonds); 
the solid line is given by constexp(— 2n/6)/n which corresponds to the solution of eq. (2). 

Figure 4. Convergence of the ground state energy in the tight-binding model of a 
finite one-dimensional lattice with A as a hopping parameter, left panels, and for a shifted 
harmonic oscillator with the hamiltonian H = a^a + X{a + a''), right panels. The upper parts 
show the energy deviation A£'„ = Eo{n) — £'0(00) as a function of the truncated dimension 
n (solid hues for A = 1 and A = 2); dotted lines show the predicted analjd;ical convergence of 
the models, Xn"^ /v? (left) and A^"/n! (right). The lower parts characterize the convergence 
rate A„ ^ A by plotting A„ = i^EnV? j-K^, left, and An = (A£'„n!)^/^"', right. 
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